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CHARACTERIZATION OF BETA DISTRIBUTION ON SYMMETRIC CONES 


BARTOSZ KOLODZIEJEK 


Abstract. In the paper we generalize the following characterization of beta distribution to the sym¬ 
metric cone setting: let X and Y be independent, non-degenerate random variables with values in (0,1), 
then U = 1 — XY and V = are independent if and only if there exist positive numbers pi, i = 1, 2, 3, 

such that X and Y follow beta distributions with parameters (pi -I- P3, P2 ) and (p 3 , pi), respectively. 


1. Introduction 

In the paper we generalize the following characterization of beta distribution to random matrices and, 
more generally, to random variables valued in the symmetric cone: let X and Y he independent, non¬ 
degenerate random variables with values in (0,1), then U = 1 — XY and V = are independent 

if and only if there exist positive numbers pi, i = 1,2,3, such that X and Y follow beta distributions 
with parameters (pi +^ 3 ,^ 2 ) and {p 3 ,pi), respectively. This univariate result was proved in [16] under 
additional assumptions that X and Y have densities, which are strictly positive on (0,1) and are log- 
locally integrable. Regularity assumption on densities was removed in the work of nnj. It turns out that 
the existence of densities assumption is redundant, what was shown in |15| . 

Here we are interested in a generalization of density versions of the beta characterization, when random 
variables are valued in the cone of r x r positive definite symmetric real matrices. Define the analogue 
of (0,1) interval in D+: Dy = {x G D+: J — x G D+}, where I is the identity matrix. Beta distribution 
on symmetric cone f2_|_ with parameters {p, q) for p,q > dim 12+ /r — 1 is defined by its density 

H(p,g)(dx) = -—2--(detxf-'i'”^/’'det(J-x)«-'i'“^+/’'/^,^(x)dx, x G D+, 

[p, q) 

where is the normalizing constant. For any x G 12+ there exists unique y G 12+ such that = x. 

Matrix y is denoted by y = x^/^. We will show that if X and Y are independent random variables valued 
in X>+, having continuous densities, which are strictly positive on X>+, then U = I — • Y ■ X^!"^ and 

V = • {I — X) ■ are independent if and only if there exist numbers pi > dimf2+/r — 1, 

i = 1,2,3, such that X and Y follow matrix-variate beta distribution with parameters (pi +P 3 ,P 2 ) and 
(p 3 ,pi), respectively. 

Actually, we will consider much more general form of transformation of random variables, which is 
defined through, so-called, multiplication algorithm. A multiplication algorithm is a mapping w: 12+ i-A 
GL{r, K) such that w(x) -w^ (x) = x for any x G 12+, where GL{r, R) is the group of invertible rxr matri¬ 
ces and ry^(x) is the transpose of w{x.). Multiplication algorithms (actually their inverses called division 
algorithms) were introduced by |13] alongside the characterization of Wishart probability distribution 
(see also [5] for generalization to symmetric cone setting). The two basic examples of multiplication 
algorithms are wi(x) = x^/^ (x^/^ being the unique positive definite symmetric square root of x) and 
W 2 (x) = 2x, where t^ is the lower triangular matrix from the Cholesky decomposition of x = fx • ij. 
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We will consider the independence of C/ = I — w{X)-Y■'w~^{X) and V = {w{U))~^ ■ {I — X) ■ , 

where w and w are two multiplication algorithms satisfying additionally some natural conditions. It turns 
out that, depending on the choice of multiplication algorithms, the characterized distribution may not 
be the beta distribution (see Theorem [6]). For example, when w = w = W 2 the condition of independence 
of U and V characterizes wider family of distributions called beta-Riesz, which include beta distribution 
as a special case. 

As in the famous Lukacs-Olkin-Rubin Theorem (see [14] for 17+ case and |3] for all symmetric cones) 
the assumption of invariance under the group of automorphisms of distributions of X and Y is considered. 
The distribution of X is said to be invariant under the group of automorphisms ii O ■ X ■ = X for 

any orthogonal matrix O. This approach leads to a characterization of beta distribution regardless of the 
choice of multiplication algorithms (see Theorem [8|). 

We cannot give the explicit formula for densities for any multiplication algorithms. In general case, the 
densities are given in terms of, so-called, zc-logarithmic Cauchy functions, that is, functions that satisfy 
the following functional equation 

/(x) -b f{w{I) ■ y • w'^{I)) = f{w{x) ■ y • w^{x)), (x, y) e f7+. 

The form of ze-logarithmic Cauchy functions without any regularity assumptions for two basic examples 
of multiplication algorithms were recently considered in [8]. Later on we will write fD(x) for the the linear 
operator acting on 17+ such that rD(x)y = w(x) • y • w^{x.). rD(x) will also be termed a multiplication 
algorithm. 

Analogous characterization of Wishart distribution, when densities of respective random variables are 
given in terms of ui-logarithmic functions is given in [Tj. Unfortunately, we cannot answer the question 
whether there exists multiplication algorithm resulting in characterizing other distribution than beta or 
beta-Riesz. Moreover, the removal of the assumption of the existence of densities remains a challenge. 

The idea of the proof is analogous to that of m- The independence condition gives us the functional 
equation for densities, which is then solved. As was observed in m, in univariate case, the independence 
condition leads to the generalized fundamental equation of information, that is 

F{x)+G[^)=H{y)+K(^^), 

where (x, y) S Dq = {{x,y) G (0,1)^: x + y € (0,1)} and F,G,H,K: (0,1) —>■ R are unknown functions. 
Our proof will heavily rely on the solution to the generalization of this equation to the cone 17+, which 
was given in |^. 

Similar characterization of beta distribution for random matrices was proved under numerous addi¬ 
tional assumptions in [^. The characterization of 2 x 2 matrix-variate beta distribution was also given 
by 0, but the characterization condition was of a different nature. 

All above considerations can be generalized to the symmetric cones, of which 17+ is the prime example. 
The paper is organized as follows. In the next section we give necessary introduction to the theory of 
symmetric cones. Next, in Section 3 we define beta and beta-Riesz probability distributions on symmetric 
cones. Main theorems are stated and proved in Section 4. Section 5 is devoted to the analysis of the 
problem, when X and Y have distributions invariant under the group of automorphisms. 

2. Preliminaries 

In this section we recall basic facts of the theory of symmetric cones, which are needed in the paper. 
For further details we refer to jl). 

A Euclidean Jordan algebra is a Euclidean space E (endowed with scalar product denoted (x,y)) 
equipped with a bilinear mapping (product) 

E X E 9 (x, y) i-A xy e E 


CHARACTERIZATION OF BETA DISTRIBUTION 


3 


and a neutral element e in E such that for all x, y, z in E: 

(i) xy = yx, 

(ii) x(xV) = x2(xy), 

(hi) xe = X, 

(iv) (x,yz) = (xy,z). 

For X G E let L(x): E —>■ E be linear map defined by 

L(x)y = xy, 

and define 

P(x) = 2L2 (x) -L(x2 ) . 

Let End(E) denote the space of endomorphisms of E. The map P: E i-A End(E) is called the quadratic 
representation of E. 

An element x is said to be invertible if there exists an element y in E such that L(x)y = e. Then y is 
called the inverse of x and is denoted by y = x“^. Note that the inverse of x is unique. It can be shown 
that X is invertible if and only if P(x) is invertible and in this case (P(x))~^ = P (x“^). 

A Euclidean Jordan algebra E is said to be simple if it is not a Cartesian product of two Euclidean 
Jordan algebras of positive dimensions. Up to linear isomorphism there are only five kinds of Euclidean 
simple Jordan algebras. Let K denote either the real numbers K, the complex ones C, quaternions HI or 
the octonions O, and write Sr(K) for the space of r x r Hermitian matrices valued in K, endowed with 
the Euclidean structure (x, y) = Trace (x • y) and with the Jordan product 

(1) xy = i(x-y + yx), 

where x • y denotes the ordinary product of matrices and y is the conjugate of y. Then S'r(K.), r > 1, 
S'r(C), r > 2, S'r(HI), r >2, and the exceptional £' 3 ( 0 ) are the first four kinds of Euclidean simple Jordan 
algebras. Note that in this case if K 7 ^ O, then 

(2) P(y)x = y-x-y. 

The fifth kind is the Euclidean space n > 2, with Jordan product 

(3) (a;o, • ■ ■, Xn) (yo, ■ • ■, J/n) = XiPi, xoyi + yoxi ,..., xopn + yox^ . 

To each Euclidean simple Jordan algebra one can attach the set of Jordan squares 

U = {x^: X G E} . 

The interior U of U is a symmetric cone. Moreover Q. is irreducible, i.e. it is not the Cartesian product 
of two convex cones. One can prove that an open convex cone is symmetric and irreducible if and only if 
it is the cone il of some Euclidean simple Jordan algebra. Each simple Jordan algebra corresponds to a 
symmetric cone, hence there exist up to linear isomorphism also only five kinds of symmetric cones. The 
cone corresponding to the Euclidean Jordan algebra equipped with Jordan product JS]) is called the 
Lorentz cone. 

We denote by G(E) the subgroup of the linear group GL{E) of linear automorphisms which preserves 
n, and we denote by G the connected component of G(E) containing the identity. Recall that if E = S'r(K.) 
and GL{r, M) is the group of invertible r xr matrices, elements of G(E) are the maps g: E —>■ E such that 
there exists a G GL(r, M) with 

g(x) = a-x-a^. 

We define K = G Cl 0(E), where 0(E) is the orthogonal group of E. It can be shown that 

(4) K = {k G G: ke = e}. 
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A multiplication algorithm is a map ^ G: x rD(x) such that tTi(x)e = x for all x G This 
concept is consistent with, so-called, division algorithm g, which was introduced by [13] and |3], that is 
a mapping 9 x i—g(x) G G such that g(x)x = e for any x € fl. If to is a multiplication algorithm 
then g = lxi“^ is a division algorithm and vice versa, if g is a division algorithm then to = g~^ is a 
multiplication algorithm. 

By in Proposition IIL4.3], for any g in the group G, 

det(gx) = (Det 5 )’’/‘^™^detx, 

where Det denotes the determinant in the space of endomorphisms on Inserting a multiplication 
algorithm g = rD(y), y € fl, and x = e we obtain 

(5) Det(lu(y)) = (dety)'^™'^/’’ 


and hence 

( 6 ) det(rD(y)x) = dety detx 

for any x, y S fl. 

One of two important examples of multiplication algorithms is the map rDi(x) = P(x^/^). The 
remaining part of this section is to give the necessary background for the definition of the second basic 
example of multiplication algorithm, the one connected with Cholesky decomposition. 

We will now introduce a very useful decomposition in E, called spectral decomposition. An element 
c S E is said to be a idempotent if cc = c 7 ^ 0. Idempotents a and b are orthogonal if ab = 0. Idempotent 
c is primitive if c is not a sum of two non-null idempotents. A complete system of primitive orthogonal 
idempotents is a set {ci,..., c,.} such that 

r 

Ci = e and CiCj = SijCi for I < * < j < r. 

i=l 

The size r of such system is a constant called the rank of E. Any element x of a Euclidean simple Jordan 
algebra can be written as x = some complete {ci,..., c,.} system of primitive orthogonal 

idempotents. The real numbers A^, i = I,..., r are the eigenvalues of x. One can then define determinant 
of X by det x = HLi 

If c is a primitive idempotent of E, the only possible eigenvalues of L(c) are 0, ^ and I. We denote by 
E(c,0), E(c, i) and E(c, I) the corresponding eigenspaces. The decomposition 

E = E(c, 0) 0 E(c, i) © E(c, I) 

is called the Peirce decomposition of E with respect to c. Note that P(c) is the orthogonal projection of 
E onto E(c, I). 

Fix a complete system of orthogonal idempotents (ci)^_^. Then for any i,j G {I,..., r} we write 

Eii = E(ci, I) = Rci, 

Ey = E n E (^Cj, iii^ j. 

It can be proved (see (U Theorem IV.2.1]) that 

E = ^ Ey 

i<j 
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and 

■ Ey c Eu + Ejj, 

Eij ■ Ejk C Ejfc, Hi ^ k, 

Eij -Eki = {0}, if {i,j}n{k,l} = 0. 

The dimension of Ey is, for any i ^ j, a. constant d called the Peirce constant. When E is S'r(K), if 
{ei,..., Cr} is an orthonormal basis of K’', then E^i = EcieJ and Ey = K (e^ej + for i < j and d 

is equal to dim\sK. 

For 1 < fc < r let Pfe be the orthogonal projection onto E*^^) = E(ci + .. .+Cfe, 1), det*-^^ the determinant 
in the subalgebra E^^\ and, for x G fl, Afe(x) = det^^^(Pfe(x)). Then Afc is called the principal minor of 
order k with respect to the Jordan frame Note that Ar(x) = detx. For s = (si, ... ,Sr) G K'’ 

and X G n, we write 

A,(x) = Ai(x)*i-*^A 2 (x)«^-'*^ ... A,(x)*g 

As is called a generalized power function. If x = X][=i then As(x) = ... a®'’. For s G K’’ 

and A G K we will write 


’ + A — (si + A,.... 


A). 


We will now introduce some basic facts about triangular group. For x and y in fl, let xDy denote the 
endomorphism of E defined by 

xDy = L(xy) + L(x)L(y) - L(y)L(x). 

If c is an idempotent and z G E(c, A) we define the Frobenius transformation Tc{z) in G by 

rc(z) = exp( 2 znc). 

Given a Jordan frame the subgroup of G, 


r= { rei(z(i))...rc,._i(z 


(r-l)'. 


E 


OiiC-i 1 ! Oil ^ 


G 0 E 


'jk 


k=j+l 


is called the triangular group corresponding to the Jordan frame For any x in fl there exists 

a unique t^ in T such that x = t^e, that is, there exist (see [U Theorem VI.3.5]) elements z^-^A g 
©L Ejk, i < j < r — 1 and positive numbers ai,... ,ar such that 


X = Tei(z(l))...Tc,_^(; 


,(r 1 )^ . 

Mapping lxi 2 : fl —> P, x M- 102 (x) = tx is the second important example of a multiplication algorithm. 


For E = 5'r.(R) we have 17 = 17+. Let us define for 1 < i, j < r matrix /iy = {'Jki)ick i<r such that 
7 y = 1 and all other entries are equal 0. Then for Jordan frame where Ck = p^kk, fc = 1,..., r, we 

have Zjk = {pjk + Pkj) & Ejk and ||zjfe|p = 2, \ < j,k < r, j ^ k. If z^*) G ©J^^+i Ey, z = I,..., r - I, 
then there exists = (a^+i,..., ar) G such that z^®) = X]J=i+iThen the Frobenius 


transformation reads below 

rc,(zh))x = J+(ah)) • X • 
where Pi(ah)) is so called Frobenius matrix: 

r 

p;(ah)) = /+ ^ ctjPji, 
j=i+l 
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ie. bellow jth one of identity matrix there is a vector a^'^\ particularly 

/I 0 0 ••• 0\ 


7^2 = 


1 

Q !3 


Vo ar 0 1/ 

It can be shown ([H Proposition VI.3.10]) that for each t G T, x e fl and s G W, 
(7) As(tx) = As(te)As(x). 

This property actually characterizes function A^ - see Theorem O 


3. Probability distributions 


The beta-Riesz distribution on symmetric cones with parameters {s,t) G M’’ x R’' for Si > {i — l)d/2, 
ti > {i — l)d/2, i = I,... ,r, (d is the Peirce constant) is defined by its density 

^77.(s, t)(dx) — - yr Qy^(x) f2/r x)/ 2 ^(x)dx, X G H, 

where I? = {x G 11: e — x G 11} is an analogue of (0,1) interval on real line and 


Bn{s,t) = 


Tn{s)Tn{t) 


rn(s + 1) 

for gamma function of symmetric cone rn(s) = 11^=1 — (j — 1)|) (see [H VII.1.1.]). 

Beta distribution on symmetric cone fl is a special case of beta-Riesz distribution for si = ... = = 

p > dim fl/r — 1 and ti = ... = = q > dim 11/r — 1 with density 

1 


6(p,g)(dx) = 


Bn{jP,q) 


(det x)P-'^™^/" det(e - x)«-‘^™^/’'/^,(x) dx, 


X G n, 


where 6o(p, q) = ro(p) := rn(p, • ■ ■ ,p)- Basic properties of beta and beta-Riesz distri¬ 

butions on 11+ are given in mm and of beta distribution on H+ in [14] . For some recent advances in 
extending beta distribution the reader is referred to m- 


4. Characterization of generalized beta distribution 

Henceforth we will denote by H an irreducible symmetric cone of rank r. The densities of generalized 
beta distributions will be given in terms of Ki-logarithmic functions, that is functions /:!!—>■]& that 
satisfies the following functional equation 

(8) /(x) -1- /(tu(e)y) = /(tu(x)y), (x,y) G H^, 

where to is a multiplication algorithm. If / is ro-logarithmic, then e-^ is said to be m-multiplicative. 
Functional equation m for ixii(x) = P(x^/^) on H+ was already considered in [T] for differentiable 
functions and in m for continuous functions on real or complex Hermitian positive definite matrices of 
rank greater than 2. Without any regularity assumptions it was solved on the Lorentz cone by [17j . The 
general forms of tui- and tTi 2 —logarithmic functions without any regularity assumptions were given in [8]. 

It should be stressed that there exists infinite number of multiplication algorithms. If to is a multipli¬ 
cation algorithm, then trivial extensions are given by rD*'^^(x) = rD(x)fc, where k € K is fixed and K is 
defined by One may consider also multiplication algorithms of the form P{x.°^)t^i- 2 c,, a G M, which 
interpolates between the two main examples: tUi (which is a = 1/2) and ru 2 (which is a = 0). In general, 
any multiplication algorithm may be written in the form ru(x) = P(x^/^)fcx, where fcx G K. 








CHARACTERIZATION OF BETA DISTRIBUTION 


7 


To define the transformation of random variables we will use two multiplication algorithms, to and to. 
Let g and g be the corresponding division algorithms, that is, g = and g = Henceforth we will 

assume that to additionally satisfies the following natural conditions 

A. to is homogeneous of degree 1, that is rD(sx) = srD(x) for any s > 0 and x e H, 

B. continuity in e, that is limx-j.e fD(x) = rD(e), 

C. surjectivity of the mapping H 9 x i-A g(x)e € H, 

D. differentiability of the mapping H 9 x i -9 rD(x), 

and the same is assumed for to. 

Conditions A — C are assumed in order to use the result of [3] regarding the generalized fundamental 
equation of information on H (see Theorem [5] below) and D is assumed to ensure that the Jacobian of 
the considered transformation exists. By tUg and tUg we will denote Iti(e) and rD(e) respectively. 

We start with the direct result, where we show that if X and Y have densities of the form (|9]), then 
the transformed variables are independent. Recall that H is an irreducible symmetric cone of rank r and 
lu-multiplicative function / satisfies the following functional equation 

/(x)/(rD(e)y) = /(rD(x)y) 

for any (x, y) € H. 

Theorem 1 Assume that multiplication algorithms tu and to are differentiable (condition D) and let X 
and Y be independent random variables valued in "D with densities of the form 

fxix) = cxidetx)‘^^’^^^''i{x)h{x)f{e- x)It>{x), 
fvix) = CYh{Wex)i{e - niex)lT,{x), 

where 

• * is to- and Vo-multiplicative, 

• f is m multiplicative, 

• h is Xo-multiplicative, 

and cx, cy are normalizing constants. Then 

U=e — m{X)Y and V = g{U){e — X) 
are independent random variables. 

Proof. Define the mapping if: T>^ —>■ T>^ by formula 

i/i(x, y) = (e - tu(x)y, g(e - tu(x)y)(e - x)), 

where g = Then we have {U,V) = if^X, Y) and the inverse mapping -0“^: TY' TY is given by 

(x,y) =i/>-^(u,v) = (e-tu(u)v,g(e-fD(u)v)(e-u)) , 

where g = Hence '0 is a bijection. We will find the Jacobian of ijj~^ in two steps. Let us observe 
that = ((2 o (fi with 

v) = (e - fD(u)v, e - u) = (a, b), 

02 (a,b) = (a,g(a)b) = (x,y). 


Denote by A the Jacobian of mapping (fi, i = 1,2. We have 


da/du 

da/dv 


da/du —rD(u) 

db/du 

db/dv 


-Idn 0 


Det(rD(u)) 
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and 


Jo — 


dx/da 

dx/db 


Ido 

0 

dy/da 

dy/db 


dy/da 

fl(a) 


= Det( 0 (a)). 


Finally, by ([5]), we get 


dim 


j = J1J2 = 


det u 


^det(e — rD(u)v) 
The joint density f(u,v) of (.U, V) is given by 


( 10 ) 


f{u,v)iu,^) = /x(e-TO(u)v)/y(fl(e-tn(u)v)(e-u)) 


det u 


dim Q 


det(e — rD(u)v) 


where fx and fy denote the densities of X and Y, respectively. Inserting m into m and re¬ 
peatedly using multiplicative properties of respective functions (that is, if h is in-multiplicative, then 
/i(x)li(rDe 5 (x)y) = /i(y) for any x, y G 11 ), we obtain 

/(; 7 ,v)(u, v) = cxcv(detu)‘^™^/’'i(u)/(u)li(e - u)/i>(u) • f{wev)i{e - Wev)Ij,{v), 
what completes the proof. □ 


Remark 1 Note that if i{x) = (det a:)^’i = (det a:)^^ dimn/r h{x) = (det dimn/r 

with Pi > dimfl/r - 1, i = 1, 2,3, then {X, Y) ~ B(pi +p^,P 2 )® B{p 3 ,pi) and {U, V) ~ B{pi +p 2 ,P 3 ) ( 8 > 
B{p 2 ,Pi), regardless of the choice of to and to. 

In order to prove the harder part of the characterization we will need the following result regarding 
the solution to fundamental equation of information on symmetric cones (see [9l Theorem 3.5]). Recall 
that V — {x € il: e — xGlI} and define 

Po = {(a,b) a + heV}. 

Theorem 2 Let a,b,c,d: > K. be continuous functions that satisfy the following functional equation 

a{x) + &(g(e - x)y) = c{y) + J(g(e - y)x), {x, y) G Vq. 

If multiplication algorithms to = g~^ and tn = g“^ satisfy conditions A — C, then there exist real constants 
Ci, f = 1,..., 4, and continuous functions hi, i = 1, 2,3, where 

• hi is XV- and xv -logarithmic, 

• h2 is XV logarithmic, 

• /13 is Xv-logarithmic, 

such that for any x G D, 

a{x) = hi{e - x) -\- h 2 {x) + h^^e - x) + Ci, 
b{x) = hi{e— XVgx) + h 3 {xvex) -|- C 2 , 
c{x) = hi{e — x) Y /i 2 (e - a;) -I- /i 3 (a:) -|- C 3 , 
d{x) = hi{e— XVgx) + h 2 {xvex) + C 4 , 

and Cl Y C 2 = C 3 Y C 4 . 

We are now ready to prove the main theorem. 
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Theorem 3 (Characterization of generalized beta distributions) Let X and Y be independent random 
variables valued in T) with continuous and strictly positive densities. Let additionally ip: TY- be a 

mapping defined through 

ip{x, y) = {e- lu(a;)y,5(e - Xo{x)y){e - x )), 

where to = and ro = g“^ are multiplication algorithms satisfying conditions A — D. Lf components of 
vector {U,V) = ip{X,Y) are independent, then there exist continuous functions i, f, g, where 

• i is to- and to-multiplicative, 

• f is to multiplicative, 

• h is to-multiplicative, 

and holds. 

Proof. Let us note that, as in the proof of Theorem[Tl the joint density of {U, V) has the form (flUl) . This 
equality is satisfied almost everywhere with respect to Lebeasgue measure. According to the assumption 
that {U,V) has independent components, we have /((7 y)(u, v) = /( 7 (u)/y(v). 

Since the respective densities are continuous, (fTUl) holds true for any u,v G V. Taking logarithm of 
both sides of (fTUl) (it is permitted, since densities are strictly positive on V), we obtain 

a (rD(u)v) + b (g(e — lu(u)v)(e — u)) = c(e — u) + d(v), (u, v) G 2?^, 

where 

a(u) = log fxie - u) - logdet(e - u), 

r 

b{u) = log /y(u), 

/ \ 1 p / \ dim n 

c(u) = log fu[e - u) -log det(e - u), 

r 

d{u) = log /y(u). 


for u G V. Let us take u = e — y and v = g(e — y)x. Then y gV, because u G T). Moreover, it is clear 
that X G n. Since v GV and e — v = g(e — y)(e — x — y) G V, we have e — (x + y) G Thus, 

( 11 ) a(x) + 6 ( 0 (e - x)y) = c(y) + d(Q{e - y)x) 

for any (x, y) G Pq = {(^-jb) G V^ : a + b G V}. Theorem [U] implies that there exist continuous functions 
hi, h 2 and /13 such that 

• hi is tu- and to-logarithmic function, 

• /i 2 is fD logarithmic, 

• /13 is tu-logarithmic. 


a(x) = lii(e - x) + /i 2 (x) + li 3 (e - x) + Cl, 
&(x) = hi{e- rOeX) + /l 3 (rDeX) + C 2 , 


for real constants Ci, i = 1,2. That is, for x G 2T we have 


( 12 ) 

(13) 


/y(x) = = eC2e^3(roeX)g/ll(e-meX)^ 


what is essentially ® for i(x) = e^^x)^ _ gii 2 (x) _ g?i 3 (x)^ 


□ 
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As was mentioned earlier, the general form of ir-logarithmic functions is known in two basic examples, 
namely to = toi = P(x^/^) and to = 032 = tx G T. These forms (see Theorem |4] and [5] below) will be 
needed in the proof of the main theorem. The proofs of these results may be found in [5] . Function H is 
called generalized logarithmic, if H{ab) = H{a) + H{b) for any positive a and b. 

Theorem 4 (to i-logarithmic Cauchy functional equation) Let f: ^ R be a function such that 

f{x) -b f{y) = / (p y) ’ y^ ^ 

Then there exists a generalized logarithmic function H such that for any x G LI, 

f{x) = H{det x). 

Theorem 5 (rD 2 -logarithmic Cauchy functional equation) Let /: —>■ R 6 e a function satisfying 

fix) + fiv) = fityx) 

for any x and y in the cone of rank r, ty G T, where T is the triangular group with respect to the 
Jordan frame Then there exist generalized logarithmic functions Hi,... ,Hr such that for any 

X G LI, 

r 

fix) = '^HkiAkix)), 

k=l 

where Afc is the principal minor of order k with respect to 

Remark 2 If we impose on f in Theorem\^ some mild conditions (eg. measurability), then there exists 
s G R*" such that for any x G LI, 

fix) = logAs(a;). 

We may now give the specification of Theorem [3l when ([9]) is known explicitly. For every generalized 
multiplication to and to, the family of generalized beta measures (as defined in ([H])) contains the beta laws 
(see Remark [T]). For to = to = toi, there are no other distributions, while for to = to = 032 generalized 
beta measures consist of the beta-Riesz distributions. It is an open question whether there is a generalized 
multiplication to that leads to other probability measures in this family. 

Define 1 := (!,...,!) G R*". Recall that LI is an irreducible symmetric cone of rank r, d is its Peirce 
constant and V = {x G LI: e — x G fl}. 

Theorem 6 (Characterization of beta and beta-Riesz distributions) Let X and Y be independent random 
variables valued in T) with continuous and strictly positive densities. Let additionally ip: ^ be a 

mapping defined through 

ipix, y) = {e- Xoix)y,Q{e - w{x)y){e - x)), 

where 03 = and to = are multiplication algorithms satisfying conditions A — D. Assume that 
components of vector (U, V) = ipiX, Y) are independent. 

If 

(1) tD(a;) = t 33 (a;) = P(a;^/^), then there exist constants pi > dimfl/r — \,i= 1,2,3, such that 

X Bipi+pz,P 2 ) and F ~ 6 (p 3 ,pi), 

(2) 03 ( 3 ;) = 03 ( 3 ;) = tx, then there exist vectors Si = (sij)j^i, stj > {j — l)d/2, i = 1,2,3, j = 1,... ,r, 
such that 


A ~ B7?.(si-b S 3 , S 2 ) and Y ^ BTZ{s 3 ,Si) 
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(3) rD(a;) = and rD(a;) = tx , then there exist constants pi > dimil/r — 1, z = 1,3 and vector 

5 2 = (s 2 ,j)j=i, 32 ,j > (j — l)d/2, such that 

XBn{{pi + ps)!, S 2 ) and YB{p 3 ,pi), 

( 4 ) tTi(a;) = tx and lxi(a;) = then there exist constants pi > dimil/r — 1, z = 1,2 and vector 

5 3 = (s 3 ,j)j=i, S 3 ,j > [j — l)d/ 2 , such that 

X ~ S72.(pi 1 + S 3 ,p 21 ) and YBTZ{s 3 ,pil). 

Proof. We start with (fT^ . If rD(x) = tni(x) = P(x^/^), then by Theorem 0] we know that there exist 
constants G R such that hi{x) = Kilogdetx, z = 1,2,3. Thus X follows B{pi +P 3 ,P 2 ) distribution 
and Y follows B{p 3 ,pi) distribution, where Pi = Ki + dimfl/r > dimfl/r — 1, z = 1, 2, 3. 

If, in turn, rD(x) = rD 2 (x) = tx, then by Theorem [5] and Remark we get the existence of vectors 
ti £ R’’ such that hi(x) = log At;(x), z = 1, 2, 3. So X follows BTZ{si + S 3 , S 2 ) distribution and Y follows 
BTZ{s 3 ,si) distribution, where Si = ti + dimll/r, z = 1,2,3 are such that Si,j > {j — l)d/2, z = 1,2,3, 
j = 1,... ,r. 

Points (3) and (4) are proved analogously. □ 

5. Distributions invariant under the group of automorphisms 

In the famous Lukacs-Olkin-Rubin Theorem (see [14] for 11+ case, [3] for all irreducible symmetric 
cones and jT] Remark 4.4] for its density version), the following independence property was analyzed: 
assume X and Y are independent random variables valued in fl and V = X+Y and U = q{X-\-Y)X (here 
g = ro'^) are also independent (supplemented with some technical assumptions). If the distribution of U 
is invariant under the group K of automorphisms, that is kU = U for any k G K, then X and Y follow 
Wishart distribution with the same scale parameter, regardless of the choice of multiplication algorithm 
in = g^^. In that case U was beta distributed for any measurable division algorithm g. Similar approach 
in our case also leads to the characterization of beta distribution on fl (see Theorem [5| below). 

Function /: fl —>■ R is called RT-invariant if /(fcx) = /(x) for any k £ K and x G 11. We will need the 
following result of [9], where compared to Theorem [2 additional assumption of RT-invariance is imposed 
on unknown functions. 

Theorem 7 Let a,b,c,d, in, in be as in Theorem 0 but assume additionally that any two unknown 
functions are K-invariant. Then, there exist constants nj, j = 1,2,3 and Ci, z = 1,..., 4, such that for 
any x £T>, 

a{x) = (ki + K 3 ) logdet(e — x) K 2 logdet a;+ Ci, 
b{x) = Ki logdet(e — x) K 3 logdet a; + C 2 , 

(14) 

c(a;) = (ki + K 2 ) logdet(e — x) K 3 logdet a; + C 3 , 
d{x) = Ki logdet(e — a;) + zt 2 logdet a; + C 4 ,, 

and Cl + C 2 = C 3 + C 4 . 

Theorem 8 (Characterization of beta distribution) Let X and Y be independent random variables valued 
in T) with continuous and strictly positive densities. Assume additionally that the distributions of X and 
Y are invariant under the group K of automorphisms. Let ip: T)^ -£ T>^ be a mapping defined through 

ip{x, y) = {e- tn(a;)z/,g(e - xo{x)y){e - x )), 

where to = g~^ and in = g“^ are multiplication algorithms satisfying conditions A — D. If components of 
vector (U,V) = ip{X,Y) are independent, then there exist constants pi > dimfl/r — 1, z = 1,2,3, such 
that X B{pi -\-p 3 ,P 2 ) and Y ~ B{p 3 ,pi). 
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B. KOLODZIEJEK 


Proof. The proof begins exactly the same as in Theorem [S] we start with (HH). If distributions of X 
and Y are invariant under the group of automorphisms, then their densities are if-invariant functions, 
that is fx{kx.) = /x(x) and /y(fcx) = /y(x) for any k G K and x S From this we conclude that 
a(u) = log /x(e — u) — disill logdet(e — u) and b{u) = log /f(u) are also if-invariant, thus by Theorem 
[3 we get the assertion. □ 
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